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SYSTEMS OF HESS-APPEL'ROT TYPE 



Vladimir Dragovic^ and Borislav Gajic^ 

Abstract. We construct higher-dimensional generalizations of the classical Hess- 
AppeProt rigid body system. We give a Lax pair with a spectral parameter leading 
to an algebro-geometric integration of this new class of systems, which is closely 
related to the integration of the Lagrange bitop performed by us recently and uses 
Mumford relation for theta divisors of double unramified coverings. Based on the 
basic properties satisfied by such a class of systems related to bi-Poisson structure, 
quasi-homogeneity, and conditions on the Kowalevski exponents, we suggest an ax- 
iomatic approach leading to what we call the "class of systems of Hess-Appel'rot 
type" . 
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1. Introduction. Starting from the Kowalevski analysis 

It is well known that Kowalevski, in her celebrated 1889 paper [28], starting with 
a careful analysis of the solutions of the Euler and the Lagrange case of rigid-body 
motion, formulated a problem of describing the parameters (A, B, C, xo,yQ, zq), for 
which the Euler - Poisson equations have a general solution in a form of uniform 
functions only with moving poles as singularities. Here, / = diag(yl, B, C) repre- 
sents the inertia operator, and x = (a;o,2/0)^o) is the centre of mass of the rigid 
body. 

Then, in §1 of [28], some necessary conditions were formulated and a new case 
was discovered, now known as Kowalevski case, as a unique possible beside the cases 
of Euler and Lagrange. However, considering the situation where all momenta of 
inertia are different, Kowalevski came to a relation analogue to the following (see 
[24]): 

^A{B - C) + yo ^/B{C - A) + zo ^/C{A - B) = 0, 

and concluded that xq = uq = zq, giving the Euler case. 

But, it was Appel'rot who noticed in the beginning of 1890's, that the last 
relation admits one more case, not mentioned by Kowalevski: 

xo^A{B -C)+ ZQ^JC{A - B) = 0, 2/0 = 0, 

under the assumption A > B > C . Such systems were considered also by Hess, even 
before Appel'rot, in 1890. Such intriguing position corresponding to the overlook in 
the Kowalevski paper, made the Hess-Appel'rot systems very attractive for leading 
Russian mathematicians from the end of XIX century. After a few years, Nekrasov 
and Lyapunov managed to provide new arguments and they demonstrated that the 
Hess-Appel'rot systems didn't satisfy the condition investigated by Kowalevski, 
which means that conclusion of § 1 of [28] was correct. 

And, from that moment, the Hess-Appel'rot systems were basically left aside, 
even in modern times, when new methods of inverse problems, Lax representations, 
finite-zone integrations were applied to almost all known classical systems, until 
very recently. 

A few years ago, we constructed a Lax representation for the Hess-Appel'rot 

system (see [15]). 

Now, in this paper the first higher - dimensional generalizations of the Hess- 
Appel'rot systems are constructed. For each dimension n > 3, we give a family of 
such generalizations. Wc provide Lax representations for all new systems, general- 
izing the Lax pair from [15]. We show that the new systems are isoholomorphic. 
This class of systems was introduced and studied in [16], in connection with the 
Lagrange bitop. 

Lax matrices of isoholomorphic systems have specific distributions of zero entries. 
Therefore standard integration techniques of [17], [1] cannot be applied directly. Its 
integration requires more detailed analysis of geometry of Prym varieties and it is 
based on Mumford's relation on theta - divisors of unramified double coverings. 

In the present paper, in addition, we perform in detail the integration procedure 
in the first higher-dimensional case n = 4 of new Hess-Appel'rot type systems. 
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The X-operator, a quadratic polynomial in A of the form A^C + AM + F, in the 
case n = 4, satisfies the condition 

Li2 = L21 = L34 = L43 = 0. 

Such situation, explicitly excluded by Adler-van Moerbeke (see [1], Theorem 1) and 

implicitly by Dubrovin (see [17], Lemma 5 and Corollary) was studied for the first 
time in [16]. (A nice and natural cohomological interpretation of polynomial Lax 
equations has been studied in [26].) 

Study of the spectral curve and the Baker- Akhiezer function for the four-dimen- 
sional Hess-Appel'rot systems shows that, similarly to [16], the dynamics of the 
system is related to a Prym variety 11. It is connected to the evolution of divisors 
of certain meromorphic differentials flj. From the condition on zeroes of the Lax; 
matrix, it follows that differentials 

arc holomorphic during the whole evolution. Compatibility of this requirement with 
dynamics is based on Mumford's relation (see [35], [16]) 

n- c e, 

where 11" is a translation of a Prym variety H, and is the thcta divisor. 

The paper is organized as follows. In Section 2, the definition of the classical 
Hess-Appel'rot system is given and a few of its basic properties are listed such as 
the A pair from [15] and the Zhukovskii geometric interpretation from [48, 31]. 
A construction of four-dimensional generalizations of the Hess-Appel'rot system is 
done in Section 3. In the same Section, a Lax representation is presented and the 
spectral curve calculated. The next Section contains generalizations of the Hess- 
Appel'rot systems to all dimensions higher than 4. In the cases n > 4 not only 
invariant relations exist, which are typical for the Hess-Appel'rot systems, but also 
values of some of the first integrals have to be fixed (and equal to zero). Thus, 
the systems we construct in the case n > 4 are also certain generalizations of the 
Goryachev-Chaplygin systems (see, for example, [24] for the definition). The Lax 
pairs are given in this section as well. In Section 5, a transformation of coordinates is 
performed for the four-dimensional Hess-Appel'rot systems, based on the decompo- 
sition ,so(4) = so(3) ® so(3). In this manner, the integration of the four-dimensional 
Hcss-Appcl'rot systems reduces to integration of two coupled three-dimensional 
systems of Hess-Appel'rot type. Starting with Section 6, the algebro-geometric in- 
tegration is performed. The principal observation is the relationship between the 
Baker-Akhiezer functions of the four-dimensional Hess-Appel'rot system and the 
Lagrange bitop. Then, in Sections 7 and 8, some most important facts from the 
algebro-geometric integration of the Lagrange bitop, done in [16], are reviewed. 
Analysis of a Prym variety H is done and through the Mumford-Dalalyan theory, 
a connection between the algebro-geometric and the classical approach from Sec- 
tion 5 is explained. Differentials are defined and the holomorphicity condition 
is derived. Therefore, the whole class of such systems is called isoholomorphic 
systems. The crucial point is application of Mumford's relation on theta- divisors 
of unramified double coverings to derive formulae in theta-functions for such sys- 
tems. Their dynamics is realized on the odd part of the generalized Jacobian, 
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which is obtainc'd by ghiing of the iiifiiiitc points of the spectral curve. In Section 
9, additional equations, which differ the cases of the Lagrange bitop and higher- 
dimensional Hess-Appel'rot systems are derived. In the final part of the paper, the 
characteristic properties, common for the Hess-Appel'rot system and its higher- 
dimensional generalizations are studied. The most relevant ones are abstracted as 
the axioms of the class of systems of Hess-Appel'rot type. In this way, in 
Section 10, after analysis of relevant Poisson structures, the Hamilton perturbation 
and bi-Poisson axioms are formulated. These axioms give very simple and geo- 
metrically transparent description of the systems of Hess-Appel'rot type. Namely, 
suppose bi-Poisson structure {•, -ji + A{-, •}2 is given, with a bihamiltonian system 
with the Hamiltonian Hq corresponding to the first structure. Further, let /i, . . . , /fe 
be the commuting integrals of the system {Hq, {•, -ji), which are Casimirs for the 
second structure {•, •}2- Then, the systems of Hess-Appel'rot type are Hamiltonian 
with respect to the first structure with a Hamiltonian 

fe 

H^Ho+Y^Jibifi, 
1=1 

where Ji are constants and bi are certain functions on the phase space. The invariant 
relations are 

fi=0, 1=1,. ..,k. 

Thus, the invariant manifolds are symplectic leaves of the second structure. 

In Section 11, a Kowalevski analysis is performed. As a result, the quasi- 
homogeneity and the arithmetic axiom are formulated, providing characterization of 
Hess-Appel'rot systems in terms of arithmetic conditions on Kowalevski exponents. 
This gives strong constraints on the functions bi in the above expressions. 

In this way, the study of Hcss-Appcl'iot systems, in a sense, roaches its historical 
origins of Kowalevski, Appel'rot, Lyapunov and others, as briefly mentioned above. 

Finally, based on these axioms we study three-dimensional Hess-Appel'rot sys- 
tems and formulate conditions which determine uniquely the classical Hess-Appel'rot 
system among them. This confirms the reasonability of the chosen axioms. Clas- 
sification of higher-dimensional Hess-Appel'rot systems looks like an interesting 
problem. We hope that detailed analysis of dynamical properties of systems of 
Hess-Appel'rot type will deserve sufficient attention. 

2. Classical Hess-Appel'rot system. 

The Euler-Poisson equations of the motion of a heavy rigid body in the moving 
frame are [24]: 

M = MxO-|-rxx, 
(1) t = Txn 

fl = JM, J = diag( Ji , J2 , J3) , 

where M is the kinetic momentum vector, fl the angular velocity, J a diagonal 
matrix, the inverse of inertia operator, F a unit vector fixed in the space and x is 
the radius vector of the centre of masses. 
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It is well known ([24]) that equations (1) have three integrals of motion: 

F, = ^{M,n) + {T,x) 

(2) F2 = (M,r) 

F3 = (r,r) = i. 

Thus, for complete integrability, one integral more is necessary [24]. Let Ji < 
J2 < J3 and X = {^OjUotZq)- Hess in [27] and Appel'rot in [4] found that if the 
inertia momenta and the radius vector of the centre of masses satisfy the conditions 

2/0 = 

(3) rz r rz r 

xoy J3 - J2 + zo\j J2 - Ji = 0, 

then the surface 

(4) Fi =M^xq + M3ZQ = 

is invariant. Integration of such a system by classical techniques can be found in 
[24]. In [15], an L-A pair for the Hess- Appel'rot system is constructed: 

L{X) = [L{X),A{X)], 

X(A) = A^C + AM + r, A{X)^Xx + ^, C = ^x, 

J2 

where skew-symmetric matrices represent vectors denoted by the same letter. Also, 
basic steps in algebro-geometric integration procedure are given in [15]. 

The Zhukovskii geometric interpretation of the conditions (3) [48, 31]. Let us 
consider the ellipsoid 

M? Ml 

-J- + ^ + ^ = 1, 
J\ J2 J3 

and the plane containing the middle axis and intersecting the ellipsoid at a circle. 

Denote by I the normal to the plane, which passes through the fixed point O. Then 
the condition (3) means that the centre of masses lies on the line I. 

Having this interpretation in mind, we choose a basis of moving frame such that 
the third axis is I, the second one is directed along the middle axis of the ellipsoid, 
and the first one is chosen according to the orientation of the orthogonal frame. In 
this basis (see [13]), the particular integral (4) becomes 

Fa = M. = 0, 



matrix J obtains the form: 

J = 



Jl 





Jl3 





Jl 





Jl3 





J3 



and X = (0, 0, Zq)- This will serve us as a motivation for a definition of the four- 
dimensional Hess-Appel'rot system. 
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3. Four-dimensional Hess-Appel'rot system. 



The Euler-Poisson equations of motion of a lieavy rigid body fixed at a point 
are Hamiltonian on the Lie algebra e(3), which is the semi-direct product of Lie 
algebras and so(3). Since is isomorphic to so(3), there are two natural 
higher- dimensional generalizations of Euler-Poisson equations. One is to Lie algebra 
e(n) = X so{n), and the second one, given by Ratiu in [37], is to the semi-direct 
product so(n) x so{n). The main result of this Section is a construction of an 
analogue of the Hess-Appel'rot system on so(n) x so(n). 

Equations of a heavy n-dimensional rigid body on so{n) x so(n), introduced by 
Ratiu in [37], are: 

M=[M,n] + [T,x] 

T = [T,n], 

where M, f2, P, x G so(n), and x is a constant matrix. We will suppose that 



n = JM + MJ, 



where J is a constant symmetric matrix. First, in this section, we consider equations 

(5) in dimension four. Motivated by the Zhukovskii geometric interpretation given 
at the end of the previous section, wc start with the following definition 

Definition 1. The four- dimensional Hess-Appel'rot system is described by the 
equations (5) and satisfies the conditions: 



(6) 



fl = MJ + JM, J = 



/ Jl 





Jl3 








Jl 





J24 


Jl3 





J3 





\ 


J24 





J3 



b) 






X12 








-X12 




















X34 








-X34 






The invariant surfaces are determined in the next lemma. 

Lemma 1. For the four- dimensional Hess-Appel'rot system, the following relations 
take place: 

M12 = Jl3(Mi3Mi2 + M24M34) J24(Mi3M34 + M12M24), 
M34 = Jl3(-Mi3M34 - M12M24) + J24(-Mi3Mi2 - M24M34). 



In particular, ifMi2 = M34 = hold at the initial moment, then the same relations 
are satisfied during the evolution in time. 
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Proof follows by direct calculations from equations (5), using (6). 
Thus, in the four-dimensional Hess-Appel'rot case, there are two invariant rela- 
tions 

(7) Mi2 = 0, M34 = 0. 

Now we will give another definition of the four-dimensional Hess-Appel'rot con- 
ditions, starting from a basis where the matrix J is diagonal in. 
Let J = diag(Ji, J2, J3, ^4)- 

Definition 1'. The four- dimensional Hess-Appel'rot system is described by the 
equations (5) and satisfies the conditions: 
a) 



b) 



X = 



c) 



+ JM, 


J = 


diag(Ji, J2, J3, J4) 





X12 


Xl4\ 


-X12 





X23 1 





-X23 


X34 ' 


-Xl4 





-X34 / 


J3- 


J4 = 


J2 — Jl, 


J3- 


Jl 





where 



ti 

t2 



2(X14X34 - X12X23) 

Xi4 - xi4 + x?2 - xis ' 

2(xi4Xi2 - X23X34) 
-X?4 - Xi4 + X?2 + xls ' 



Proposition 1. There exists a bi- correspondence between sets of data from Defi- 
nition 1 and Definition 1 '. 

Proof. From J = S'^JS, where 

/ cos(p sint^ 



S = 



cosyJi siny3i 
sin if cos ip 
\ — sin cos (fix , 



and <^ = 5 arctan (^1 = 5 arctan ^^77, we have 

J = diag 



J1 + J3- A J1 + J3-A1 J1 + J3 + A J1 + J3 + A1 



Here A = ^(-^3 - ^1)2+4723, Ai = ^(-^3 - Ji)^ + 4JI4 . The first part in the 
Definition I'c) 

J3 — Ji = J2 — Jl, 
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follows from these relations. 
Prom X = S'^xS, we have 








Xl2 





Xl4 




-Xl2 





X23 





[ 





-X23 





X34 


-Xl4 





-X34 






where 

Xi2 = Xi2 cos ipcostpi+ X34 sin (p sin ipi , 
Xi4 = Xi2 cos (fi sin (fii - X34 sin ip cos (pi, 
X23 = -Xi2 sin (f cos </Ji + X34 COS sin , 
X34 = Xi2 sin ipsirnpi+ xm cos cos i^i • 

Prom the last formulae, it follows: 

(8) (xi2 sin (fi + X23 cos (p) cos v?i + (xi4 sin (p - X34 cos (p) sin = 0, 



(9) (xi2 cos (fi - X23 sin (p) sin (pi - (xm cos (p + X34 sin cos ipi = 0. 

From (8) and (9): 

. o 2(xi4X;m - X12X23) , 
13.11 Zip = To To tt— =: tl, 

Xf4 - X34 + XI2 - X23 



, o 2(xi4Xi2 - X23X34) , 
tan2yi = — ^2 — —To — —To- =: i2- 

-Xi4 - X34 + XI2 + X23 



Thus, we get 



(Ji - J3) 



2 



(J3 - Ji)^ 

1 + 4 ' 



2 _ ( J4 - J2) 



2 



From the last formulae, we come to the last part of the Definition 1' c. This finishes 
the proof. □ 

Note. 1) In the case 

■^24 7^ 0, X34 = 0, 

there is an additional relation 

X12X34 + X14X23 = 0. 

It follows from the system 

X12 sin <p + X23 cos ip = 0, 
Xi4 sin p} - X34 cos (p = 0, 



as a consequence of (8, 9). 
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2) In the case 
additional relations are 



J24 = 0, X34 = 0, 



X34 = Xi4 = 0, 

and the second relation from the Definition 1' c) can be replaced by the relation 



X12 Y J2- Ji + X23 V J3- J2 = 0. 

Notice the similarity of the last condition with the condition (3) for the classical 
three-dimensional case. (By ignoring the last coordinate one can recover the three- 
dimensional Hess-Appel'rot case.) 

Theorem 1. The four- dimensional Hess-Appel'rot system has the following Lax 
representation 



i(A) = [L(A),A(A)], 
i(A) = A^C + AM + r, A{X) = Xx + il, C = 



1 



X- 



Proof. Proof follows from 
where 



J1 + J3 









Dl3 

































/ 



Dl3 = -Xl2( J13M12 + J24M34) + X34( J13M34 + J24M12), 
D24 = -xMJlsMsi + J24M12) + X34{Jl3Mi2 + J24M34), 

using relations (7). □ 

One can calculate the spectral polynomial for the four-dimensional Hess-Appel'rot 
system: 

p(A, = det(L(A) - M • 1) = + -P(A)m^ + QiXf, 



where 



P(A) = aA^ + bX^ -\-cX'^ + dX + e 
Q{X) = fx"" + gX^ + hX" + iA + i 

fl = (^12 + ^"34, 

b = 2C12M12 + 2C34M34{= 0), 

c = + M^^ + MI3 + MI4 + M^2 + + '2C12T12 + 2C34r34, 

d = 2T12M12 + 2ri3Mi3 + 2ri4Mi4 + 2r23M23 + 2r24M24 + 2r34M34 

6 = r^2 + + r^4 + r23 + r24 -i- t\^, 

f = C12C34 

g = Ci2M3i + C3iMi2{=0), 

h = T34C12 + T12C34 + M12M34 + M23M14 - M13M24, 

i = M34T12 + M12T34 + Mi4r23 + M23ri4 - ri3M24 - r24Mi3, 

j = r34ri2 + r23ri4 — ri3r24. 



10 



VLADIMIR DRAGOVIC AND BORISLAV GAJIC 



Let us consider standard Poisson structure on scmidirect product ,so(4) x so(4). 
The functions d,e,i,j are Casimir functions (see [37]), c,h are first integrals, and 
b = 0,g = are the invariant relations. General orbits of co-adjoint action are eight- 
dimensional, thus for complete integrability one needs four independent integrals 
in involution. 



4. The n-dimensional Hess-Appel'rot systems. 



In this Section, we introduce Hess-Appel'rot systems of arbitrary dimension. 



Definition 2. The n-dimensional Hess-Appel'rot system is described by the equa- 
tions (5), and satisfies the conditions: 

a) 



(10) 



n = JM + Mj, J = 



/ Jl 





Jl3 





. 







Jl 





J24 


. 


. 


^13 





Ja 





. 


. 





J24 





Js 


. 


. 














. 


. 


V 











. 





b) 



X = 








X12 


. 






-X12 





. 


. 










. 


. 










. 


. 


v 








. 


'. 0/ 



Direct calculations from (5) using (6) give the following lemma: 



Lemma 2. For the n-dimensional Hess-Appel'rot system, the following relations 
are satisfied: 
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a) 



Mi2 = Jl3(Mi2Mi3 + M24M34 + ^ M2pM3p) + 



p=5 



J24(Mi2M24 + M13M34 - ^ MipM4p) 



p=5 



M34 = -Jl3(Mi3M34 + M24M12 + ^MipMp4)- 

p=5 



J24(Mi3Mi2 + M24M34 + M2pM3p), 

Msp = - Jl3(Mi3M3p + M2pMi2) - J24(M34M2p + M2sMip)- 

n 

Msi^ip - Q.zAMip + ^(M3fef^fep - fi3fcM4p), p > 4, 

fe=5 

Map = Ji3(-Mi4M3p + M1PM34) + J2i{Mi2Mip - M2iMip)- 

n 

Mu^sp + naiMsp + ^(M4fcf^/cp - fiikMip), p > 4, 



b) 



fc=5 



Mm = 0, A: J > 4. 



c) T/ie n-dimensional Hess-Appel'rot case has the following system of invariant 

relations 

(11) Mi2 = 0, Mip=0, l,p>3. 

By diagonalizing the matrix J, we come to another definition 

Definition 2'. The n-dimensional Hess-Appel'rot system is described by the equa- 
tions (5), and satisfies the conditions 
a) 

fl = JM + MJ, J = diag(Ji, J2, J3, J4, J4), 

b) 



X ■ 



X12 





Xl4 •• 


• '\ 


X12 


X23 


.. 


. 


-X23 





X34 •• 


. 


Xl4 


-X34 


.. 


. 








.. 




J3 — J4 










J4- 


J2 










X12X34 + X14X23 = 
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where 



_ 2(xi4X34 - X12X23) 
X?4 - Xu + X?2 - xis ' 

2(xi4Xi2 - X23X34) 



t2 ■■ = 



-Xli - X34 + X12 + X23 " 



As in the dimension four, there is an equivalence of the definitions. 

Proposition 2. There exists a hi- correspondence between sets of data from Defi- 
nition 2 and Definition 2'. 

Proof follows the steps in Proposition 1. 

Next theorem gives a Lax pair for the n-dimensional Hess-Appel'rot system. 
Theorem 2. The n-dimensional Hess-Appel'rot system has the following Lax pair 

i(A) = [L(A),A(A)], 



i(A) = A^C + AM + r, A{\) = \x + ^, C- 



1 



Jl + J3 



X- 



Proof. The statement follows from 



[C,0] + [M,x] 



X12 



J13M12 + J24M34 J24M45 

J13M34 + J24M12 J13M45 





V 



J24Af4n \ 





/ 



and relations (11). □ 

Note. Let us note that invariant relations (11) exist in a more general case, 
with matrix J given by: 



J:= 



/ Jl 





Jl3 


Jl4 


Jin \ 





Jl 


J23 


J24 ■ 


■ J2n 


Jl3 


J23 


J3 


. 


. 


Jl4 


J24 





J3 . 


. 


V Jin 


J2n 





. 


: J3/ 



But, using transformations J 1— > T* JT, where T is a block-diagonal matrix with 
2 X 2 - block A € S0{2) and (n - 2) x (n - 2)-block B € SO{n-2) on the diagonal, 
such a more general matrix J can be transformed to the case considered above. 

Let us mention again that the Goryachev-Chaplygin system is a classical case 
integrable for the fixed level of a first integral. According to Lema 2b) and 2c) in 
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the n-dimensional Hess-Appcl'rot system we also fix values of certain first integrals. 
But also, we have invariant relations which do not exist in the Goryachev-Chaplygin 
case. 



5. The decomposition so(4) = so(3) © so(3) and 
integration of the four-dimensional Hess-Appel'rot system. 

Starting from the well-known decomposition so(4) = so(3) ® so(3), let us intro- 
duce 

Ml = ^(M+ + M_) M2 = ^(M+ - M_), 

(and similarly for ^l,T,x), where M+,M_ are vectors in defined with follow- 
ing correspondence between two three-dimensional vectors and four-dimensional 
antisymmetric matrices 



(M+,M_) 



/ 

Ml 

-Ml 



V Mi 

Then, equations of the motion become 



-Ml Ml -Ml \ 
-M2 
-Mi 
J 





M| 

M2 



-Ml 


M3 



(12) 



Ml = 2(Mi X Oi + Ti X xi) 
M2 = 2(M2 X O2 + X X2) 



f 1 = 2(ri X Qi) 

r2 = 2(r2 X ^2), 



and 



1 



1 



Xi = (0, 0, --(xi2 + X34)), X2 = (0, 0, --(xi2 - X34))- 



Integrals of the motion are 

(M„M)+2 

(13) 



1 



J1 + J3 



) = hi, 

(^^,^,) = l, 
(Mi,r,) = c„ 

(Xi,Mi) =0. 



1,2, 



Connections between M and are 



f2l = ((Ji + J3)M(i)i - (Ji3 - J24)M(2)3, (Jl + J3)M(i)2, 
(Jl + J3)M(i)3 (Jl - J3)M(2)3 - ( Jl3 + J24)M(2)i), 

il2 = ((Jl + J3)M(2)i - (Ji3 + J24)M(i)3, (Ji + J3)M(2)2, 
(Jl + J3)M(2)3 + (Jl - J3)M(i)3 - (Ji3 - J24)M(i)i), 



where Mj^^^ is the j-th component of vector Mj. Using these expressions, equations 
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(12) can be rewritten in the following form: 

M(i)i =2[( Ji - J3)M(i)2M(2)3 - (Jl3 + J24)M(i)2M(2)i + r(i)2X(i)3], 
M(i)2 =2[-(Ji - J3)M(2)3M(i)i - (Ji3 - J24)M(i)3M(2)3 + 

(Jl3 + '/24)M(i)iM(2)l - r(i)iX(i)3], 
M(l)3 =2(Ji3 - J24)M(i)2M(2)3, 

(14) f (1)1 =2[r(i)2((Jl + J3)M(i)3 + (Ji - J3)M(2)3 - (Jl3 + J24)M(2)i)- 

r(l)3(>/l + ^3)M(i)2], 
r(l)2 =2[r(i)3((Ji + J3)M(i)i - (Ji3 - J24)M(2)3)- 

r(i)i((Ji + J3)M(i)3 + (Ji - J3)M(2)3 - (Jl3 + ./24)M(2)i)], 
f (1)3 =2[r(i)i(Ji + J3)M(i)2 - r(i)2((Jl + J3)M(i)i - (Ji3 - J24)M(2)3)], 

and 

M(2)l =2[(Ji - J3)M(2)2M(i)3 - (Ji3 - J24)M(2)2M(i)i + r(2)2X(2)3], 
Af(2)2 =2[-(Jl - J3)M(i)3M(2)i - (Ji3 + J24)M(2)3M(i)3 + 

(Jl3 - J24)M(2)iM(i)i - r(2)iX(2)3]! 
M(2)3 =2(Ji3 + J24)M(2)2M(i)3, 

(15) f (2)1 =2[r(2)2((Jl + J3)M(2)3 + (Jl - J3)M(i)3 - (J13 - J24)M(i)i)- 

r(2)3(^l + ^3)M(2)2], 
f (2)2 =2[r(2)3((Jl + J3)M(2)i - (Ji3 + J24)M(i)3)- 

r(2)l((Jl + J3)M(2)3 + (Jl - J3)M(i)3 - (Ji3 - J24)M(i)i)], 
r(2)3 =2[r(2)i(Jl + J3)M(2)2 - r(2)2((Jl + J3)M(2)i - (J13 + J24)M(i)3)]. 

Prom the equations (14) and (15), it follows that M(i)3 = M(2)3 = 0, giving two 
invariant relations introduced before. 

Now, we are going to proceed the integration in a classical manner. 
First, let us introduce coordinates Ki and U as follows: 

Af(j)i = sin k, M(i)2 = Ki cos k, i = l,2. 

From the sixth equation of (14), using integrals (13), we have that 

2 



r^i)3 = 4(Ji + J3)' 



(1 - r^i)3)(/ii - -^_^x(i)3r(i)3) - 4 



= ^3(r(i)3). 



Thus r(i)3 can be solved by an elliptic quadrature. Also from the energy integral 
(the first one in (13)) we have that 

^1 = - :^>^(i)3r(i)3- 



Since tanZi = using first two equations in (14), we have: 

2X(i)3Ci 



h = -2(Ji3 + J24)-ftr2 sin/2 + 
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Also from the second and third integral in (13), we have that 

Kfrl^^ - 2ciM(i)2r(i)2 + cf - Mfi)i(l - Tl^,) = 0. 
Similarly, from equations (15), we get: 

2 

(1 - r(2)3)(/l2 - J.J X(2)3r(2)3) " C2 



^2)3 = 4(Jl + J3)' 



^3(r(2)3), 



.2 - 2 



K2 = h2- -J— -J-X(2)3r(2)3, 

12 = -2(Ji3 - J24)K, sin/i + , 

^|r^2)2 - 2C2M(2)2r(2)2 + - Mf2)i(l - r^2)3) = 0- 

Prom the previous considerations, we conclude that for complete integration 

of the four-dimcnsional Hess-Appcl'rot system one need to solve a system of two 
differential equations (for h and h) of the first order and to calculate two elliptic 
integrals, associated with elliptic curves Ei and E2 defined by 

(16) Ei-. = Pi{x) = SAiX^ - 4,Bix'^ - 8AiX - 4Ci, i=l,2 

where 

Ai = (Jl + J3)X(i)3, Bi = (Ji + J^fhi, d = (Ji + Jsfic^i - hi). 

This is a typical situation for the Hcss-Appcl'rot systems that additional integra- 
tions are required (see [36, 24, 15, 13]). Now we pass to the algebro-geometric 
integration. 



6. Algebro-geometric integration 



Before analyzing spectral properties of the matrices L{X), we will change the 
coordinates in order to diagonalize the matrix C. In this new basis the matrices 
L{X) have the form Z(A) = U-'^L{X)U, where 



U 



/ 








i\/2 
2 


f \ 










V2 






2 


2 






%/2 










2 


2 


\ 


2 


iV2 
2 





/ 



After straightforward calculations, we have 
/ 



L(A) 



iA34 

ZA34 

/?3 - iBi -id?. + Bi 



-PI - iPl 

-101 - PI 

-?Ai2 





Z/33-/34 \ 

-/33 + ipi 


«Ai2 



where 



A12 = A^Ci2 + AM12 +ri2, 
A34 = X^Csi + AM34 + r34. 
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(17) 



/33 


= xz + Ays, 


X3 = 


I 

2 (ri3 + ir2s), 


Pi 


= Xi + Xy4, 


X4 = 


i (ri4 + ir24) , 


PI 


= X3 + Ays, 


2/3 = 


^ (Mis + iMss) , 


PI 


= X4 + Xy4, 


2/4 = 


^ (Mi4 + iM24) . 



Matrix L(A) is of the same form as the Lax matrix for the Lagrange bitop [15, 
16]. It is a quadratic polynomial in the spectral parameter A with matrix coef- 
ficients. General theories describing the isospectral deformations for polynomials 
with matrix coefficients were developed by Dubrovin [17, 18] in the middle of 70's 
and by Adler and van Moerbeke [1] a few years later. Dubrovin's approach was 
based on the Baker- Akhiezcr function. Both approaches were applied in rigid body 
problems (see [32, 1] respectively). 

But, as it was shown in [16], none of these two theories can be directly applied in 
cases like this. Necessary modifications were suggested in [16], where a procedure 
of algebro-geometric integration was presented. It is based on some nontrivial facts 
from the theory of Prym varieties, such as the Mumford relation on theta divisors 
of unramified double coverings and the Mumford-Dalalyan theory (see [16, 35, 34, 
14, 40, 5]). 

Here, we are going to follow closely the procedure from [16], with necessary 
changes, calculations and comments. 

As usual, we start with the spectral curve 

r : det (l(A) - /i • l) = 0. 

We have 
(18) 

r : + (A?2 + A34 + 4/33/3s* + ^^PiPl) + [A12AS4 + 2i{p;p4 - p^P^f = 0. 
There is an involution 

a : (A,/L<) (A, -/^) 

of the curve F, which corresponds to the skew-symmetricity of the matrix L{X). 
Denote the factor-curve by Fi = F/cr. 

Lemma 3. The curve Fi is a smooth hyperelliptic curve of the genus g{Ti) = 3. 
The arithmetic genus of the curve F is (?a(F) = 9. 

Proof. The curve 

Fi : + P{X)u + [Q(A)]2 = 0, 
is hyperelliptic, and its equation in the canonical form is: 

(19) ui = ra! _ [Q(A)]^ 

where wi =u + P{X)/2. Since - [Q(A)]2 is a polynomial of the degree 8, the 

genus of the curve Fi is .g(Fi) = 3. F is a double covering of Fi and the ramification 
divisor is of degree 8. According to the Riemann-Hurwitz formula, 5a (F) =9. □ 
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Lemma 4. In generic case the spectral curve T has four ordinary double points 
Si, i = 1, ... ,4:. The genus of its normalization T is five. □ 

Proof. From the equations 

dp{X, m) ^ Q dp{X, /^) ^ Q 
dX ' dn ' 

where p(A, fi) = det (^L{X) — fi ■ ij — fi^ + fj,^P{X) + [(5(A)]^, the double points are 

Sk = (Afe,0), k = 1, ... ,4, where Afe are zeroes of (5(A). Thus, g{T) = (?a(r) — 4 = 
5. □ 

Lemma 5. Singular points Si of the curve T are fixed by a. The involution a 

exchanges the two branches of T at Si. 

Proof. Fixed points of the a are defined with /U = 0, thus Si are fixed. Since their 
projections on Fi are smooth points, a exchanges the branches of F, which are 
given by the equation: 

, _ -P(A) + ^P^{X) - 4Q2(A) 
^ ~ 2 

□ 

We start with the well-known eigen-problem 



(20) + A{X)j = 0, L(A)^fe = ^iki^k, 

where are eigenvectors with eigenvalues [ik ■ Then ipk (t, A) form a 4 x 4 matrix 
with components ili\{t,X). Denote by its inverse matrix. Let us introduce 

g]{t,{\^ik))=Mt,X)■^%X) 

(there is no summation on k) or, in other words g{t) = tpkit) ® fit)''. 

Matrix g is of rank 1, and we have dip/dt = —Atp, d(p/dt = (pA, dg/dt = 
[g,A\. We can consider vector-functions 4'k{t,X) = (^'^(f. A), V'fc(t, A))'^ as one 
function tP{t, (A, fi)) = {^p^{t, (A, ^)), ...,ip^{t, (A, /x)))^ on F defined by (A, /ife)) 
= ipKtjX). Similarly, we define (fi{t, {X, /j,)). Relations for divisors of zeroes and 
poles of functions i (fi in the afHne part of F are: 



(21) {g^^ = dj{t)+d'{t)-Dr-Di, 

where dj{t) is divisor of zeroes of tpj, divisor d'(t) is divisor of zeroes of t^', Dr is the 
ramification divisor over A plane (sec [17]), D'^ is some subdivisor of Dg divisor of 
singular points defined by (21). One can easily calculate deg Dr = lQ,degDs = 8. 

Matrix elements gj{t, {X,^k)) are meromorphic functions on F. We need their 
asymptotics in neighbourhoods of points Pk, which cover the point A = oo. Let V-'fe 
be the eigenvector of the matrix L{X) normalized in Pk by the condition i/:^. = 1, 
and let (^*^ be the inverse matrix for t/'^. We will also use another decomposition 
of matrix elements of g: gj = tpl^j = '4'kV'j- It is an immediate consequence of 
proportionality of the vectors ipk and ipk {^p'^ and (p'^). 
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Lemma 6. 

a) Matrix g has the following representation 



^ dp{X,ij)/dfjL 

where ai = L, a2 = P • 1 + i^, as = PL + L^. 
b) For the Lax matrix L and Aj such that Q{Xi) = 0, it holds as = 0. 

The proof of the Lemma follows from [17] and straightforward calculation. From 
the part (a) one can see that g could have poles in singular points of the spectral 
curve. But, from (b) we have 

Corollary 1. The matrix g has no poles in singular points of the curve F. 

So, from now on, taking Corollary 1 into account, wo will consider all functions 
in this section as functions on the normalization T of the curve F. 

Since the functions -0^ and (p'^ are meromorphic in neighbourhoods of points Pk, 
their asymptotics can be calculated by expanding ipk as a power series in A~^ in a 
neighbourhood of the point A = oo around the vector Ck, where e^ = 51. We get 



(22) 




where matrices C, M and F are defined by L(A) = X^C + AM + F. Comparing the 
same powers of A, from (22) we get 

{ui)i = 0, {vi)i = 0, {wi)i = 



(23) 



{Ui)j 
{Wi)j 



M. 









1 




~ ^33 




1 


Ci 




= Mi, 


, d. 



3 7^ i 



MjkMki 



MiiMji ~ 
+ 1 1 



Cii — Ckk Cii — c. 



33 



^Mjk{vi)k + jfc(ui)fe - bi{vi)j - di{ui)j 



y-v Mtfc^i +fit, hi = ^Mik{vi)k +^f'jkiui)k 



Cii Ckk 



kjLi ^1% ^KK k^i k^i 

So, the matrix ip = {^|} in a neighbourhood of A = oo has the form: 

1 

A3. 



(24) 



u 
X 



V w 

A2 + A3 



O 



Denote by dj and rf' the following divisors: 

dl=di+P2, d2=d2 + Pl, d3 = d3 + P4, d4 = d4 + P3, 

d^ = d^+P2, d^ = d^ + Pi, d^ = d^ + Pi, d^ = d'^ + P3. 
Analyzing the behavior of matrix g around points Pk, as in [16], we get 
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Proposition 3. 

a) Divisors of matrix elements of g are 

(25) (5}) = +dj-Dr + 2 (Pi +P2 + P3 + Pi) -Pi- Pj 

b) Divisors di,d^ are of the same degree 

degdi = deg d^ =5. 

Let us denote by $(t, A) the normalized fundamental solution of 
|+i(A)U(t,A)=0, $(r)=l. 



.9t 

Then, if we introduce the Baker- Akhiezer functions 

(26) ^\t, T, (A, Mfc)) = E ^sit, A)/i^(r, (A, /z,)) 

s 

where are eigen-vectors of L(A) normalized by the condition J2s ^^{^^ Mfc)) = 
1, it follows that: 

Proposition 4. Functions tp^ satisfy the following properties 

a) In the affine part of T, the function tp^ has 4 time dependent zeroes which belong 
to the divisor d^{t) defined by formula (21), and 8 time independent poles, i.e. 

(i,\t,T,{\,iJ.k)))^=d\t)-V, degP = 8. 

b) In points Pk, functions tp^ have essential singularities as follows: 

r{t, T, (A, /x)) = exp [-{t - T)Rk - iFu)] a'{t, r, (A, /x)) 
where Rk and Fk are: 

Ri = i (^) , R2 = -Ri, R3 = I (^) , Ri = -R3, 
Fi = (^J^ Qaidt^ , F2 = -Fi, F3 = (^J^ ni2dt^ ,Fi = -F3 

and are holomorphic in a neighbourhood of Pk, 

d'(T, T, (A, 11)) = h\r, (A, a\t, r, Pk) = 6^ + 4 + 0{z'), 

with 
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Proof repeats the demonstration of Proposition 5 in [16]. 

Let us denote by t/i^^ the Baker- Akhiezer function for the Lagrange bitop from 
[16] with analytical properties as in Proposition 4 a) above and with asymptotics 
given by: 

at points Pfe, functions tjj\^g have essential singularities as follows: 
i^LBit'.'r, (A,/u)) =exp [-(i - T)i?fc] ^'^^(t, r, (A,/i)), 
where Rk are given with 

R\=i {^~^ > -^2 = —Ri, R3 = i (^^) ' ~ 
and a]^Q are holomorphic in a neighbourhood of Pk, 

dis(r,r,(A,M)) = /iis(r,(A,M)), dls(i, r, Pfe) = + + 0(^2). 

Prom the Proposition 4 and from Proposition 5 of [16], we have 

Corollary 2. A relationship between the data of generalized Hess-Appel'rot prob- 
lem and the Lagrange bitop are given by: 

V'ha := V''' = exp(iPfe)t/>is, A; = 1, . . . , 4; 

b) 

(29) v^^HA ■■= = exp(i(Pfc + Fj))v^, k,j = l,...,4. 

{vj we will also denote as Vjj^^.) 



7. A Prym variety 

Let us recall that (t) is divisor defined in (21). 
Lemma 7. On the Jacobian Jac(r) the following relation takes place : 

A{d'{t) + ad^{t)) = Aid^ir) + ad^r)) 

where A is the Abel map from the curve F to Jac(r), and a is involution on T. 

The proof is the same as the one of the corresponding Lemma in [16]. 

From the previous Lemma, we see that vectors A{d^{t)) belong to some transla- 
tion of a Prym variety 11 = Prym(r|Fi). More details concerning Prym varieties 
one can find in [41, 40, 21, 9, 34, 35, 5, 8]. A natural question arises to compare 
two-dimensional tori 11 and Ei x E2, where elliptic curves Ei are defined by (16). 

Together with the curve Fi, one can consider curves Ci and C2 defined by the 
equations 

(30) C^:v' = ^ + Q{X), C2:v' = ^-Q{X). 
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Lemma 8. Curves Ei defined by (16) are Jacobians of curves Ci given by (30). 
Proof. Follows by a straightforward calculation. 

□ 

Since the curve Fi is hyper-elliptic, in a study of the Prym variety 11 the 
Mumford-Dalalyan theory can be applied (see [14, 34, 40]). Thus, the previous 
Lemma allows us to use the following Theorem from [16]. 

Theorem 3. 

a) The Prymian 11 is isomorphic to the product of curves : 

n = Jac{Cx) X Jac{C2). 

b) The curve F is the desingularization o/Fi Xpi C2 and Ci Xpi Fi. 

c) The canonical polarization divisor S 0/ IT satisfies 

S = £1 X 62 + 61 X E2, 

where Oj is the theta- divisor of Ei. 

Theorem 3 explains the connection between the curves Ei , E^ and the Prym vari- 
ety n. Further analysis of properties of Prym varieties necessary for understanding 
the dynamics of the Lagrange bitop will be done in the next section. 

8. Isoholomorphisity condition, Mumford's relation and 
solutions for v'^j^^ 

We saw that integration of the four-dimensional Hess-Appel'rot system is par- 
tially reduced to solutions of the Lagrange bitop. Now, we arc goint to give the 
explicit formulae for the Baker-Akhiezer function for the Lagrange bitop, obtained 
in [16]. According to Proposition 4, the Baker-Akhiezer function satisfies usual 
conditions of normalized (n=)4-point function on a curve of genus .9 5 with the 
divisor V of degree degV = g + n— 1 = S, see [19, 18]. By the general theory, it 
should determine the whole dynamics uniquely. 

Let us consider the differentials fJ^ = gijdX, i,j = 1,...,4. In the case of 
general position it was proved by Dubrovin that 51* is a meromorphic differential 
having poles at Pi and Pj, with residues Vj and —vf respectively. But here we have 

Proposition 5. [16] Differentials ^i, ^4, ^3 are holomorphic during the 
whole evolution. 



The proof is based on the fact that from the conditions = -^1 = -^4 = 
it follows that 



(31) vl =vl = vl=vt = 0. 

(For more details see [16]). We can say that the condition L\ = L\ — L\ = L\ = {) 
implies isoholomorphicity. Let us recall the general formulae for v from [18] : 

i mAiP,)-A{Pj)+tU + zo) 
^ ' '~ XAtU + zo)e{Pi,Pj) ' 
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where U = ^x'^'^^U^^^ is a certain linear combination of 6-periods C/*^'^ of differen- 
tials of the s 
scalars, and 



tials of the second kind Op.'*, which have a pole of order two at P,; A, are nonzero 



<PuPj)--- 



(Here v is an arbitrary odd non-dcgcnerate characteristic.) Thus, from (32) we get: 
Holomorphicity of some of the differentials O] implies that the theta divisor of 
the spectral curve contains some tori. 

In a case when the spectral curve is a double unramified covering 

7r:f^ri; 

with ^(ri) = g, g(r) = 2g — 1, as we have here (assuming that F is the normal- 
ization of the spectral curve F), it is really satisfied that the theta divisor contains 
a torus, see [35]. Let us denote by 11" the set 

n- = {l e Pic^^-'^t\NmL = Kn, h^{L) is odd} , 

where K-p^ is the canonical class of the curve Fi and Nm : PicT — + PicTi is the 
norm map, see [35, 40] for details. For us, it is crucial that H" is a translate of the 
Prym variety 11 and that Mumford's relation ([35]) holds: 

(33) n- c Op. 
Let us denote 

(34) U = i{X34U^^^ - X34f/^'^ + Xi2C/^'^ - X12U^^^), 

where is the vector of 6-periods of the differential of the second kind Op,-* , which 
is normalized by the condition that a-periods are zero. We suppose here that the 
cycles a, b on the curve F and a, b on Fi are chosen to correspond to the involution 
a and the projection tt, see [5, 40]: 

7r(ao) = ao; 7r(6o) = 26o, o-{dk) = ak+2, k = l,2. 

The basis of normalized holomorphic differentials [uq, . . . , Ms] on F and [vq, vi, U2] 
on Fi are chosen such that 

Tr*{vo)=Uo, Tr*{vi) =Ui + a{ui) =Ui + Ui+2, i = 1,2. 

Now wc have 
Theorem 4. [16] 

a) // vector zq in (32) corresponds to the translation of the Prym variety 11 to Il~ , 
and vector U is defined by (34) then conditions (31) are satisfied. 
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b) The explicit formula for zq is 

(35) Zo = 1:{t00,T01,T02,T01,T02), TQi = j ui, i = 0,l,2. 



Formulae for scalars Aj from (32) will be given later in this section. 

The evolution on the Jacobian of the spectral curve Jac(r) gives a possibility to 
reconstruct the cvohition of Lax matrix L{X) only up to a conjugation by diago- 
nal matrices. To overcome this problem, we are going to consider, together with 
Dubrovin, a generalized Jacobian, obtained by gluing together the infinite points. 
Those points are Pi, P2, -P3, Pi and the corresponding Jacobian will be denoted by 

Jac(f|{Pi,P2,P3,P4}). 

The generalized Jacobian can be understood as a set of classes of relative equiva- 
lence among the divisors on F of a certain degree. Two divisors of the same degree 
Di and D2 are called equivalent relative to points Pi, P2, P3, Pa if there exists a 
mcromorphic function / on F such that (/) = Di — D2 and /(Pi) = /(P2) = 
fiP:,)=f{P,). 

The generalized Abel map is defined with 

i(P) = (A(P),Ai(P),...,A4(P)), Ai(P)=exp T Op,Q„, i = 1, 4, 

JPo 

where A is the standard Abel map. Here ^p^Qo denotes the normalized differential 
of the third kind, with poles at Pi and at an arbitrary fixed point Qq. 
Then the generalized Abel theorem (see [21]) can be formulated as 

Lemma 9 (the generalized Abel theorem). Divisors Di and D2 are equivalent 
relative Pi , P2 , P3 , P4 if and only if there exist integer-valued vectors N, M such that 

A{Di) = A{D2) + 2nN + BM, 

Xj{Di) = cXj{D2)exp{M,A{D2)), j = 1,...,4 

where c is some constant and B is the period m,atrix of the curve F. 

A generalized Jacobi inverse problem can be formulated as a question of finding, 
for given z, points Qi, • • • , Qs such that 

8 4 

1 2 

\j = cexp V / ^P^Qa + Kj, j = 1, ...4, 
s=l -^^0 

where K is the Ricmann constant and constants kj depend on F, points Pi, P2, P3, P4 
and the choice of local parameters around them. 

We will denote by Qs the points which belong to the divisor V from Proposition 
4, and by E the prime form from [21]. Then we have 
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■J' 



Proposition 6. Scalars Xj from formula (32) are given with 

8 fQs 

Xj = A° exp ^ la;^^) , A° = c exp ^ / Qp. q„ + 

where x = (a;(^\ . . . , a;^^)) = i(x34, -X34, X12, -X12) and 

7/ = ^^n£;(P„P)|p=p^.. 

(k^^ is a local parameter around Pj.) 

To give formulae for the Baker- Akhiezcr function, we need some notations. Let 

a^i,) = exp[^E7^-«]^(^^J|WT^, 
where ^ 

JPo 

and 7™ is defined by the expansion 

fiW = -km + 7;;: + P -> P„. 



Denote 



{x,F)-a(,x)e^V[ ^j^^l^x UpJ e{A{P) ^ A{P,) - z,) 

Finally we come to 

Proposition 7. [16] The Baker- Akhiezer function is given by 

^o e{A{P-P,)-zo) 

V^-(f,P)=^-(f,P)-^^-^^^P^, i = l,...,4, 

Z^fe=l e(P,Pk) 

where zq is given by (35). 

9. The restrictively integrable part 
equations for the functions Fi,i = 1, . . . , 4 

Let us denote 

01—P1 + P3, h ■■= Fi - F3; 

and also 

Ni := Mi4 - M23, iV;5 := Af24 - M13, 
7V2 := -M24 - Mi3, iV4 := Mi4 + M23. 

From (29) we have 

'Pi ■= arg{vlHA) =4>i+ ai{t), ip2 := arg(w]^^) = (f>2 + a2{t), 

where ai{t) = arg(v^]^g) and a2{t) = arg{vlj^g) are known function of time. Let 
us denote Uj = tant^j. 

Basic relationships among those quantities are given in the next proposition. 
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Proposition 8. The following relations take place 
a) 

b) 

(37) 7Vi = -2|i;iis(Cn -033)1 sin ^1, 7V4 = 2|viiB(Cn - C;4)| cos ^2; 

c) 

(38) ^1 =iVl(J24 + J13), ^2 = -A^4(J24 - J13). 

Proof, a) follows from the formulae for L (28) and (29). Part b) also uses Corollary 
2 b, Proposition 4 b. Note that from the condition = JM+MJ and the invariant 
relations we have: 

= M14J24 + M32J13; 
^34 ^ M32J24 + M14J13. 

From the last relation and the definition of functions from the Proposition 4, c) 
follows. □ 

Using formulae (37), (38) we get 

^1 = -2(J24 + Ji3)\{Cii - C33)t'3LBl sin((/)i + ai(t)) 

(t>2 = -2(J24 - ^13)|(C'll - C44)vIi^b\ C0s((/>2 + a2{t)) 

10. Restrictive integrability in an abstract Poisson algebra settings. 
Bihamiltonian structures for the Lagrange bitop 
and n-dimensional Lagrange top 

From the analysis given in this paper, it follows that the Hess-Appel'rot system 
and its generalizations can be understood as natural examples of the following, 
more abstract situation. 

Suppose a Poisson manifold (M^", {•, •}) is given, together with k + 1 functions 
H,fu...,fk& C°°{M), such that 
(Al) 

k 

{H,fi} = Y,aijfj, aijGC^{M), i,j = l,...,k; 

(A2) 

{,/.,./,} = 0, i,j = l,...,k. 

The Hamiltonian system (M", H) will be called restrictively integrable, if it sat- 
isfies the axioms (A1-A2). 

A more general case can be obtained by replacing condition (A2) with 
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(A2') 

k 
1 = 1 

In this case, the algebra of invariant relations is a noncommutative Lie algebra. 

Starting from the Hamiltonian system (M, Hq) with k integrals in involution 
fi, . . . , fk, choosing functions bj G C°°{M), j = 1, . . . , fc, one comes to a restric- 
tivcly integrable system: 
(HP) Hamiltonian perturbation 

The system (M, H) where 

k 

will be called a Hamiltonian perturbation. It satisfies (Al) with 

O'ij = {bj , fi} , i,j = l,...,k. 

Natural question is the converse one: when a restrictively integrable system is of 
the form (HP)? 
Denote 

From the Jacobi identity, and involutivity of functions fi we get compatibility 

conditions. 

Proposition 9. // a restrictively integrable system which satisfies the axioms (Al- 
A2) is of the form (HP), then 

^il ~ j; ^ = 1 • • • ) k. 

If in Proposition 9 we replace axiom A2 with A2' then c^j should satisiy 

^il ~ ^li~^ 53 ^il^jm- 
m 



A three-dimensional Lagrange top is defined by the Hamiltonian: 

according to the standard Poisson structure 

{Mi,Mj}i = -CijkMk, {Mi,Tj}i = -eijkTk, {Ti,rj} = 

on the Lie algebra e(3). It is also well-known that three-dimensional Lagrange 
top is Hamiltonian in another Poisson structure, compatible with first one. This 
structure is defined by: 

{ri,r,}2 = -cijkTk, {Mi,M2}2 = 1, 
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and the corresponding Hamiltonian is: 

HL = {a- 1)M3 ( i(M2 + + Fa ) + MiFi + MaFa + MsFg 

where Ji = 1, /a = a, = 1. 

Casimir functions in the second structure are + + F3 and M3. 

Let us observe that the Hamihonian for the three-dimensional Hess-Appel'rot 
case is a quadratic deformation of Hamiltonian Hl of the Lagrange top: 

Hha = Hl + J13M1M3. 

The function M3, which gives the invariant relation for the Hess-Appel'rot case, is 
a Casimir function of the second Poisson structure. 

Having this observation in mind, next we arc going to prove that the Lagrange 
bitop and the n-dimensional Lagrange top are also bihamiltonian systems. 

The standard Poisson structure on the semi-direct product so(4) x so(4) is: 

{Mij,Mjk}i = -Mik, {My,Fjfc}i = -Fifc, {Fi,-,Ffc;}i =0. 

Now let us introduce a new Poisson structure as follows: 

{Fy, Fjfc}2 = — Fjfc, {Mjj, Ffei}2 = 0, 

(39) {Mi3, M23}2 = -X12, {Mi4, M24}2 = -Xl2, 

{Mi3, Mi4}2 = -X34, {M23, M24}2 = -X34 

Casimir functions in this structure are M12, M34, Ff2+Fi3-|-Ff4-|-F|3-|-F|4-|-F§4, 

and Fi2r34 + F23ri4 - Fi3r24. 

Proposition 10. The Poisson structure (39) is compatible with the standard one. 

Proof. Two Poisson structures, defined with antisymmetric matrices A and B, are 
compatible if their Shouten bracket, defined by: 

cyclic fori, 

vanishes (see [24]). Proof follows by direct calculation. □ 

In the metric = JM + MJ, where J = diag(Ji, Ji, J3, J3), Hamiltonian func- 
tion of the Lagrange bitop in the standard Poisson structure is: 

Hlb = i(2JiMf2 + (Ji + J3)Mi^ + (Ji + J3)Mi^+ 

(Jl + J3)M|3 + (Ji + J3)M|4 + 2J3MI4) + Xl2ri2 + X34r34- 



Let us assume that Ji = a, J3 = 1 — a. 
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Proposition 11. The Lagrange bitop defined in the first Poisson structure by the 
Hamiltonian Hlb is a Hamiltonian system in the second Poisson structure (39) 
with the Hamiltonian: 



Hlb = 

(2a-l)(xi2Mi2 + X34M34 ) ( Mi^3 + Mi^4 + MI3 + MI4 

X?2 - xIa 



Xl2^l2 + X34r34 



^ (l-2a)(xi2M34 + X34M,2) (^^^^^^ _ ^^^^^^ ^ ^^^^^^ ^ ^^^^^^^^ 
X12 X34 

Mi2ri2 + Mi3ri3 + Mi4ri4 + M23r23 + M24r24 + M34r34. 

The situation with four-dimensional Hess-Appel'rot case is similar to the three- 
dimensional case: the Hamiltonian for the four-dimensional Hess-Appel'rot system 
in the first structure is again a quadratic deformation of Hlb' 

Hha = Hlb + Ji3(-Afi2M23 + M14M34) + J24(Mi2Mi4 - M23M34) 

Functions M12 and M34, giving invariant relations for the four-dimensional Hess- 
Appel'rot system, are also Casimir functions for the second Poisson structure (39). 

Putting X34 = 0, and assuming xi2 = 1 in (39) and in expression for Hlb, we 
get the bihamiltonian structure for the four-dimensional Lagrange top introduced 
by Ratiu in [37]. 

In general, in arbitrary dimension n, the standard Poisson structure on so{n) x 
so{n) is given by: 

{Mij,Mjk}i = -Mik, {Mij,Tjk}i = -Tik, {Tij,Tki}i=(i, i,j,k = l,...,n. 

In the metric O = JM + MJ, the n-dimensional Lagrange top is defined with a 
Hamiltonian 



xi2r 



Hl = U 2 JiM^^ + ( Ji + J3) J2^M^P + ^2p) + 2 J3 X] j 

y p=3 3<p<q<n J 

The number of nontrivial integrals of the motion is — [^] . Casimir functions 

are given by (see [3 7]) 

(40) tr{T^''), ir(Mr2^+i). 

We use the notation Ji = a, J3 = 1 — a, X12 = 1- Let us introduce a new 
Poisson structure: 

(41) {Tij,Tjk}2 = -Tik, {Mij,Mki}2 = {Mu,M2i}2 = -1, / = 3,...,n. 
The dimension of a symplectic leaf in this structure is 

(n-2)(n-3) 



+ 4(n-2), 



2 

hence there are " -^n+s _|_ Casimir functions: 

Mi2,Mpg,Tr(r2'=), 2<p<q<n, k = 1, . 
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Proposition 12. The n- dimensional Lagrange top is a Hamiltonian system in the 
Poisson structure (41), compatible with the standard one. Its Hamiltonian is: 

= (2a - l)Mi2 p{Ml + Ml) + T,^ + 

(l-2a) ^ Mpg(MiqM2p-M2gMip + rpg)+ ^ M^^Tpq 

S<p<q<n l^P<g^^ 

Similarly as in dimension 3 and 4, Hamiltonian for the Hess-AppeProt system 
in arbitrary dimension n is a quadratic deformation of the Hamiltonian for the 
n-dimensional Lagrange top: 

n 

HhA =Hl + ^(JlsMifeMsfe + J24M2fcM4fc), 

fe=l 

and functions M12, Mpq,p, q> 3, which give the invariant relations (11), are Casimir 
functions for the Poisson structure (41). 

We can summarize the discussion of this section by saying that constructed 
Hess-Appel'rot systems satisfy the following. 

(BP) (bi-Poisson condition) There exist a pair of compatible Poisson struc- 
tures, such that the system is Hamiltonian with respect to the first structure, having 
the Hamiltonian of the form (HP), such that fi are Casimir functions with respect 
to the second structure. 

The invariant relations define symplectic leaves with respect to the second struc- 
ture, and the system is Hamiltonian with respect the first one. 

11. Back to Kowalevski properties 

As we mentioned in the introduction, from the first years of its history, the Hcss- 
Appel'rot systems were closely related to Kowalevski's analysis. Investigating the 
systems constructed in the first part of this paper, we have noticed that they are cer- 
tain perturbations of the form (HP)of integrable systems, n-dimensional Lagrange 
tops and the Lagrange bitop. We also observed that perturbing functions fi, which 
give the invariant relations, arc Casimir functions of the second Poisson structure 
and the integrable systems are bihamiltonian corresponding to that structure. Up 
to now there is no restriction on the choice of perturbing functions bi in (HP). In 
order to define more precisely a class of systems which has the same typical dynam- 
ical and analytical properties as the classical and the n-dimensional Hess-Appel'rot 
case we need to study them in more details. Finally, after that we will be able to 
extract the basic ones leading to the constrains on the functions bi. The correctness 
of our choice is illustrated by the Theorem 5 in Section 12. Using the axioms one 
can easily construct large number of new examples of systems of Hess-Appel'rot 
type (beside semidirect product so{n) x so{n) in a study of generalized rigid body 
systems, one can consider, for example, semidirect product of i?" and so{n)). 

To get the right choice of axioms, we have to turn back to the Kowalevski 
analysis. First, we are going to introduce some general notions, see [29]. 
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Suppose a system of ODEs of the form 

(42) Zi = fi{zi,...,Zn), i = l,...,n, 

is given and there exist positive integers gi, i = 1, . . . ,n, such that 
fi{a^^zi,...,a^"'Zn) =a^*~^'^fi{zi,...,Zn), i = l,...,n. 

Then the system (42) is quasi-homogeneous and numbers gi are exponents of quasi- 
homogeneity. Then, for any complex solution C = (Ci, . . . , C„) of the system of 
algebraic equations: 

(43) -giCi = fi{Ci,...,Cn), i = l,...,n, 
one can define the Kowalevski matrix K = K{C) = [-Rr*(C)]: 

^(C) = ||(C) + 5,<5i. 

Eigen-vahies of the Kowalevski matrix arc called the Kowalevski exponents. This 
terminology was introduced in [45]. In last twenty years, heuristic and theoretical 
methods in application of Kowalevski matrix and Kowalevski exponents in study of 
intcgrability and nonintcgrability have been actively developing, see for example [2, 
3, 46, 47, 25, 29]. But the notion of Kowalevski matrix and Kowalevski exponents 
were introduced by Kowalevski herself in [28]. The criterion she used [28, p. 183, 
1. 15-22] to detect a system which is now known as the Kowalevski top, can be 
formulated in Yoshida terminology as: 

Kowalevski condition (Kc). The 6x6 Kowalevski matrix should have five 
different positive integer Kowalevski exponents. 

Now we return to the study of Hess-Appel'rot systems. The systems we have 
constructed are quasi-homogeneous. Exponents of each M variable are g = 1, 
and for any F they are equal to two. We are going now to calculate Kowalevski 
exponents for the Hess-Appel'rot systems. 

Three-dimensional Hess-Appel'rot case. Let us denote (Mi, M2, M3, Fi, 
Ts) by (^1, . . . , Zq). Then the Euler-Poisson equations take the form (42) with 

/i = (-^3 - ^1)22-23 + Ji3-zi2;2 + 2:5; 

/2 = (^3 — Jl)ziZi + ^13(23 — z1) — Z4,; 

h = -^132^22:3; 

/4 = ^3^5^3 - ■^1^2^6 + ■^13^1 2:5; 

fb = —-hzzZi + JiZiZq + Jisiz^z^ — Z1Z4); 
/e = J3Z2Z4, — J1Z1Z5 — J13Z3Z5; 

and g,; = 1, z = 1, 2, 3 and gi = 2, z = 4,5,6. The invariant relation corresponds 
to the constraint C3 = 0. So, we are looking for solutions (ci, C2, 0, C4, C5, cq) of the 
system of the form (43). One can easily get C4 = — Ji3cf + C2, C5 = — ci(l + J13C2), 
C6 = — (cf + C2)/2. Then, for c\ ^ 0, we get four possible solutions for (ci,C2) 
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divided into two pairs: (±i/Ji3, — I/J13) and (±2i/Ji3, — 2/J13). The Kowalevski 
exponents are 

(-1,-2,2,4,3,3), (-1,1,3,2,2,2), 

respectively. 

Thus, it can easily be seen that classical Hess-Appel'rot system doesn't satisfy 
exactly the Kowalevski condition (Kc), although it is quite close to. 

Four-dimensional Hess-Appel'rot systems. In the four-dimensional case, 
there are 12 variables (M12, M13, M14, M23, M24, M34, Tis, ri3, ri4, r23, r24, r34). 
Denote them by (zi, . . . , Z12) and 

fi = Mzi,...,zi2), i = l,...,12, 

the corresponding right sides of the Euler-Poisson equations. Exponents of quasi- 
homogeneity are 

gi = l, i = l,...,6, gi = '2, z = 7, ...,12. 

We are looking for solutions (di, . . . , ^12) of a system of the form (43). The invariant 
relations correspond to constraints 

di = de = 0. 

From the relations on /2, /a, /4, /s one can express (ds, ^9; dio, dn) as functions of 
(d2,...,d5). Then, from the equation on /7,/i2 one gets d7,di2 as functions of 
(d2, . . . , ds). Consider two possible cases separately. 

Example 1. First consider the case where one of perturbing constants is equal 
to zero, say J13 = 0. The solution of the last system of four equations on four 
unknowns (d2, . . . , ds), leads to the solution: 

d2 = -d5(= dg = dio), d3 = d4 = ^-1 - d^(= dn = -dg), (dy = di2 = 1), 

with arbitrary d2. Computing the Kowalevski matrix, we get finally the Kowalevski 
exponents 

(0, -1, 3, 4, 2, 1, 2,1, 2 + 2^J|4(l+di), 2 - 2yfjlJl + 4), 
1 + 2^J|4(1 + 4), 1 - 2^Jl{l + 4)). 

By analyzing correspondent eigen- vectors, we see that eight of them are tangent 
to the symplectic leaf, and other four are transversal to the leaf. Nonintegral 
Kowalevski exponents 

(2 + 2^J|4(l+d2), 2 - 2^J|4(l + di), 1 + 2^J24(l+d2), 1 - 2^J24(l + di)) 
correspond to a half of tangential eigen- vectors. 
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Example 2. Now, suppose that both Ji.j and J24 arc nonzero. To simphfy the 
computations, assume X12 = 1, X34 = 2. If d2 = 0, then there are three cases of 
nontrivial solutions of the system (43): 

1) da = 0, ch = 0, d4 = 0, 4 = ±§; 

2) d2 = 0, ds = Tj, di = ±^/4, 4 = 0; 

3) d2 = 0, (i3=±i, di = ±l, d5=0. 

Let us calculate Kowalevski exponents in the last case with ^3 = j/4, ^4 = i/4. 
We get first 

, (Ji + J3) , * ^ n ^ n ^ * ^ (-/i + J3) 
«7 = ,ds = ~2'^ ^ ^ ' " ^ 2 ^ 48 ' 

and then the Kowalevski exponents 

(0, -1, 3, ,^^+(Jii_M,^_(Jii:_M, 

,^{Ji3_M,2-^-:^l^^,2, 1, 2, 1). 

Suppose now that d2 is arbitrary. Then there are two sets of solutions of (43) of 
the form 

4) 

ct2 = s, ds = Ti- J , ^ , d4 = ±z-!^ , d5 = s; 

■Ji + Js Ji + J3 



5) 



"2 = s, ds = ±1 , , , di = ±1 , d5 = -s. 

Ji + J3 Ji + J3 

In the case 4) we get further 



, 1 , Vl + ,s2(Ji+J3)2 

dr = — J , T ' ^8 = r~~j ' = 

Ji + J3 -h + ./3 

, , . Vl + ,s2(./i+ J3)2 1 

dw = -s, du = , di2 - 



J1 + J3 Ji + J 3 

The Kowalevski exponents are 

(0, -1, 3, 4, 2, 2, 1, 1, 1 + A 1- A 2 + A, 2- A) 

where 

A = 2(Ji3 + J2iWl + s'^{Ji + J3y- 

Five-dimensional Hess-Appel'rot systems. In this case, there are 20 vari- 
ables (M12, Mi3, Mi4, Mi5, M23, M24, M25, Af34, M35, M45, ri2, ri3, ri4, ri5, r23, r24, 
r25, r34, Fss, r45). As before, we denote them by (2:1, ... , 220)- Denoting also by 

fi = fi{zi,.--,Z2o)., i = l,...,20. 
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the corresponding right sides of the Euler-Poisson equations. Exponents of quasi- 
homogeneity are 

gi = l, i = l,...,10, gi = 2, z = ll,...,20. 

We are looking for solutions (rfi, . . . , ^20) of a system of the form (43). The invariant 
relations correspond to constraints 

di=d8 = dg = diQ = 0. 

Prom the relations on /2 — /? one can express (rfi2, dn) as functions of (^2, dr). 
Then, from /n one gets dn as a function of {d2, ■ ■ ■ jd-/), and after that, from 
(/18, /i9, /20) one gets dig, dig, ^20 as functions of (^2, . • . , ^7). The final step is 
solution of the system (/12, . . . , /17) with the unknowns (^2, • . • , dj). 

Example 3. We describe nonzero solutions of (43) under the assumption ^5 = 
dr = 0, X12 = 1. There are eight sets of solutions: 
1) 

2i 

d2 =0, ds = 0, di = 0, (is = ±- 



2) 
3) 
4) 
5) 
6) 
7) 
8) 



"Ji + J3' 



2i 

±-— — -, ds = 0, d4 = 0, de = 0; 
Ji + Ja 



d2 = ±— —, ds = 0, ^4 = 0, de = ± 



J1 + J3 J1 + J3 



2i 2 
d2 = ± — , ds = 0, d4 = 0, de = - — ; 
•J13 •Jis 



d2 = — — , da = 0, d4 = ± — , de 

^13 J13 



4 4i 
d2 = — ^— , da = 0, d4 = ±-— , de = ±- 



-^13 -^13 -^1 + -^3 



1 i 2i 

d2 = -j—, ds = 0, d4 = ±-— , de = ±- 



■^13 -^13 Ji + J3 



d2 = -^,ds = -^, d, = ± ^f-^^'^\ d, = Q. 

^13 ^"^24 ^>^13>^24 



We calculate the Kowalevski exponents in three representative cases: 1), 4) and 

8). 

In the case 1) only nonzero di, z > 15 arc die — 2/(Ji + J3), dig — — 2i/(Ji + J3), 
and (under the assumption J3 = 1, J13 = 10) the Kowalevski exponents are 

(l + A, 1-A A -A B,l-B, 1, -1, -1, -1, 3, 3, 3, 2, 2, 2, 2, 4, 4, 4); 
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where A = 2 ^^"° and B = 2i^^. 

In the case 2) dn = 2/(Ji + J3) and dis = 2i/(Ji + J3) are the only nonzero 
di, i> 11 and the Kowalevski exponents are 

(-1, -1, -1, 1, 2, 2, 2, 2, 3, 3, 3, 4, 4, 4, 1 + A, 1 - A, A, -A, B, 1 - B), 

where A = B = ^ and Ji = l,Ji3 = 10, J24 = 4. 

In the case 4) we get that all dj, i> 11 are zero except 

^2 2i 

"15 = , "17 = — • 
Jl3 Jl3 

The Kowalevski exponents are 

(-2, 0, 0, 4, -1, 4, -1, 4, -1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3). 
In the case 5) the Kowalevski exponents are 

(-1, 0, 0, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3). 
In the case 8) we have the following nonzero di, i > 11: 

ai6 - j^, di7 - T j^, rfi8 - T . . . • 

4Ji3 4,724 4J13J24 

Corresponding Kowalevski exponents are 
(0, -1, I 4, ^, \/2, -s/2, 1 + ^/2, 1 - x/2, 1, i 3, ^, 1, i 3, ^, 2, 2, 2). 

Six-dimensional Hess-Appel'rot systems. In this case, there are 30 vari- 
ables (M12, . . . , M56, ri2, . . . , Tse). Denoting them by (zi, . . . , Z30), by 

fi = fi{zi,...,Zso), i = l,...,30, 

the corresponding right sides of the Euler-Poisson equations and exponents of quasi- 
homogeneity by 

gi = l, i = l,...,15, 5i = 2, i = 16, ...,30, 

we search to solutions (d\, . . . , dzo) of a system of the form (43). The invariant 
relations correspond to constraints 

di = diQ = dn = di2 = dis = du = di5 = 0. 

The solution of the system follows the same lines as in the five-dimensional case. 
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Example 4. Under the following assumptions 

Jl =1, J3 = 3, Xl2 = 1, d6 = ■■■ = dg - 

we get six sets of solutions of the system (43): 
1) 

d2 = 0, da = -, ^4 = 0, ^5 = 0; 

2) 



d2 = ds = ^rr^, d4 = 0, ds = 



13 ^'^24. 



2Ji3 2J24 V 



24 



3) 



4) 



5) 



6) 



2 ^ n ^ . /s2j23 + 4 ^ 

- — , da = 0, d4, = id j^- , d5 = s; 



— , d3 = 0, CI4 = — -^^^ — , ds = s; 



^2 = da = — , d4 = 1/ . ds = 0; 

^Jia ZJ24 V •^13^24 



3 7 / 72 _ Q 72 _ 4 2 72 t2 

. _ J _ J _ , / -^13 ^"'24 ■^24-^13 J _ o. 

'^'--Ws''^'-W4''^'-y M 



24 

where s is an arbitrary parameter. 

In the case 1) the only nonzero d are die = 5 and d2i = 5. The Kowalevski 
exponents are 

(-1, -1, -1, -1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 

4, 4, 4, 4, 1 + ^, 1 - A, A, -A, 1-iB, 1 - iS, zB, iS), 

where A = ^^^^"'^^ , B = J13/2, Ji = 1,J3 = 3, J24 = 4. 

Thus, using into account properties of Kowalevski exponents of algebraically- 

intcgrablc Hamiltonian systems, wc can conclude that for the systems we have 
constructed, functions bi in the perturbation formula (HP) should satisfy two con- 
ditions: 

(QH) (quasi-homogeneity) The obtained system of Hamiltonian equations 
has to be quasi-homogeneous. 

In such Kowalevski matrix exists and we come to the last condition. 

Suppose the invariant relations correspond to equations Zi = 0, . . . , Zk =0. 

Denote by p number of Casimirs: n = p + 2m, where 2m is the dimension of a 
general symplectic leaf. 
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(ArA) (Arithmetic axiom) For any nonzero solution C = (0, 0, Ck+i, Cn) 
of the system (43), the Kowalevski matrix K(C) has n — p eigen-vectors tangent to 
the symplectic leaf and p transversal to it. Half of the Kowalevski exponents which 
correspond to tangential eigen-vectors and all of transversal ones are rational num- 
bers. Irrational numbers among the second half of tangential Kowalevski exponents 
are divided into pairs such that the differences are integrally dependent. 

12. Description of three-dimensional systems of Hess-Appel'rot type 

Now we would like to derive conditions which determine classical Hess-Appel'rot 
system among three-dimensional systems of Hess-Appel'rot type. More precisely, 
suppose a system is given by the Hamiltonian 

(44) Hi=Ho + JbMa, 

where Hq is the Hamiltonian of the Lagrange top corresponding to the first Poisson 
structure, M3 is its integral and a Casimir for the second structure, J is a nonzero 
constant and fo is a function, such that the axioms of the systems of Hess-Appel'rot 
are satisfied. 

Analyzing the system (43) we come to the first, simple but very important, 
properties of such functions b. Denote by / the value of the function b at the point 
(ci, . . . , ce) of nonzero solution of (43). 

Lemma 10. For a nonzero solution (ci, cq) of {43) and the value f — b{ci, cq) 
it holds 

a) cl + cl=0 or f = 0; 

b) if f o,i^d ci = ±ic2 then f = Ti/J or f = ^2i/J. 

The Lemma follows by straightforward calculation. It gives a possibility to 

reduce the analysis of functions b to analysis of their germs. By a germ of a 
function b, we mean (/, /i, • • • , /e), where fi is the value of bz^ calculated at points 
of solutions of system (43). 

Lemma 10 leads to important simplifications in a study of Kowalevski matrices 
and their characteristic polynomials. We will denote by Ki,...,K4 Kowalevski 
matrices evaluated on germs, where 



Ki 


(/ = 


-i/J,ci = 


-- iC2), 


K2 


(/ = 


2i/J,ci = 


-iC2), 


K3 


(/ = 


i/J, ci = - 


-iC2), 


Ki 


(/ = 


-2i/J,ci 


= iC2), 



and by Pchj, i = 1, . . . , 4 the corresponding characteristic polynomials, Pchi(w) = 
dv\.{K^ - wid). 

Proposition 13. 

a) Characteristic polynomials Pch^ have integer-valued coefficients. 

b) These coefficients are J -independent. 

c) The characteristic polynomial Pchi is of the form 



(45) 



Pchi(u;) =w^ + Air,uf' + Auw"^ -|- • • • -|- ^10, 
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where 

(46) Ai5 = (-9 - 2 J/1C2), Aio = mhJc2{-ihJc2 + /1/2 J'c^ - iJ^clfl - i). 

The proof of Proposition 13 follows from Arithmetic axiom, Lemma 10 and 

straightforward calculations. 

Proposition 14. For i = 1, . . . , 4 the following relation holds 

Pchi(-l) = 0. 

Proposition 14 is a well-l?;nown property of Kowalevski matrices for authonomous 
systems, see [29]. For i = l, using the notation 

X := JC2/2, Y := JC2/1, 

we get the following 

Corollary. For the first germ, the following relation holds 

(47) {2XY - SiF + X - 3i) (-y + X - 2) = 0. 

In the same notations, from Proposition 13 we get 
Lemma 11. IfY^O, then 

(48) {X-i){Y-i{X + i)) = 0. 

By systematical analysis of equations (45-48) finally we come to the following 

Theorem 5. The only non-zero polynomials b which give systems of Hess-Appel'rot 
type by relation (44) are of the form 

b{zi, .. .,zq) = zi+ kz3. 

All systems of Hess-Appel'rot type of the form (44) are the classical Hess-Appel'rot 

systems. 

Example. One of possible solutions of the system (46, 47) is X = i, Y = —3. 

It leads to the function 

(49) b{zi,...,ze) = -3zi-i-iz2. 
Corresponding characteristic polynomial of the Kowalevski matrix Ki is 

Pchi(w) = w{w - l){w - 2){w - 3)(w + l){w + 2). 
However, the characteristic polynomial of the Kowalevski matrix K3 is of the form 
Pch3(u;) = {w- l){w - 2){w - 3)(«; + l){2w'^ -2w-\- 9). 
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Thus, the function b given by (49) only partially satisfies the Arithmetic axiom. 
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